NA5A-CR- 192277 


C ^ /p Aj 7 
/*S-G 7- £ 
/* 8/0 / 


p.r? 

Final Report of Phase-II Research 

on 


Applications of Active Adaptive Noise 
Control to Jet Engines 


Prepared by: 


Rahmat Shoureshi, Professor 
Larry Brackney, PhD Candidate 
Ray W. Herrick Labs 
School of Mechanical Engineering 
Purdue University 
(317)-494-5639 


Submitted to: 


M. M. Bright 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, OH 


Grant No.: NAG3-1272 

Grant Period: March 1 , 1992 through February 28, 1993 


o 

eg 


h* 



o 

rg 

ro 

iH 

l 


00 

rn 

o 


a* 

c 

(H 

z 

D 

o 



G3/07 


o o . - 

K W C 
to ms 
Z JZ 
O D © 
*-« QC *-• 3 
K K -D 

< Z * v. 

U O +> 3 

m o k- a. 

*J 0 w 

a oj a 
a on a) 

< M QC fO 

o a 

z f o 


^ ^ I— < 

Mil C 
N > — • 
rg u. jq 
rg i— a> 
^ a ^ it 


i 

QC 

U 

I 

< 

< 

z 


< LU 

O Z 00 

< t m N 
O 

uu Z t 
> UJ 
*~f tv 

>- K o 
uaio^ 

< ^ 


a 


EXECUTIVE SUMMARY: 


During phase-Il research on the application of active noise control to jet engines, the 
development of multiple-input/multiple-output (MIMO) active adaptive noise control 
algorithms and acoustic/controls models for turbofan engines were considered. Specific 
goals for this research phase included: 

♦ Implementation of a MIMO adaptive minimum variance active noise controller 

• Turbofan engine model development 

A minimum variance control law for adaptive active noise control has been developed, 
simulated, and implemented for single-input/single-output (SISO) systems [1-41- Since 
acoustic systems tend to be distributed, multiple sensors and actuators are more 
appropriate. As such, the SISO minimum variance controller was extended to the MIMO 
case. Simulation and experimental results are presented. 

A state-space model of a simplified gas turbine engine is developed using the bond 
graph technique [5,6]. The model retains important system behavior, yet is of low enough 
order to be useful for controller design. Expansion of the model to include multiple stages 
and spools is also discussed. 


MIMO MINIMUM VARIANCE CONTROLLER DEVELOPMENT: 


Before discussing the adaptive MIMO minimum variance controller, it is important to 
recall some important results from [1-4]. First, an AutoRegressive Moving Average with 
exogenous input (ARMAX) model structure is assumed for the system: 


y(t) = 


B(z) 

A(z) 


z^ u(t) + 


C(z) 

A(z) 


w(t) 


where: A(z), B(z), and C( z) are polynomials in the unit delay operator, z; d is the time 
delay between the input, u(t), and the output, y(t); and w(t) is a white noise disturbance to 
the system. 

The coefficients of the polynomials A(z) and B(z) may be identified using a Recursive 
Least Mean Squares (RLMS) algorithm. This simple iterative method consists of the 
following relations: 

a,(k+ 1 ) = a,(k) + u a (y(t) - y(t)) y(t-d-i) 
bj(k+l ) = bj(k) + % (y(t) - y(t)) u(t-d-i) 


where: the a, and bj’s are the estimated coefficients of the polynomials; u a and Ub are the 
“learning” gains: and Y is the estimated system output, given previous inputs, outputs and 
polynomial estimates. 


Finally, a SISO minimum variance control law is given as: 


where: p>0. 


u(t) = - 


1+ fX 


y(t) 


The previous results are easily extended to the MIMO case. For multiple inputs and 
outputs, the ARMAX structure now becomes: 
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In block diagram form, this relationship for a two-input/two-output system is shown as 
figure one. 



Figure one: Block Diagram of a Two-Input/Two-Output ARMAX 
Active Noise Control System 



The multivariate RLMS algorithm is now expressed as a set of matrix iterations: 
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Lastly, the MIMO minimum variance control law is: 

Hi 


u(t) = 


0 

1 + Pi 

0 

0 0 


y(t) 


1 + Hn 


A digital computer simulation has been developed to evaluate the MIMO adaptive 
minimum variance control structure. The simulated plant has the following ARMAX 
model: 


r - .2 z + .02 -,-2 

.1 z 2 - .02 z + .002 ,-2 1 
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The initial model estimate is given as: 


y(t) = 


z~ - .5 z + .5 z -2 
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z : - .5 z + .5 
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As a test of the RLMS algorithm, an initial simulation is run with w(t)=0, and uj(t) and 
U 2 (t) equal to the following input: 

1.4 [l.3 sin (4.5 n t) + .4 sin (3 • 4.5 ;i t) + .25 sin (5 • 4.5 n t)] t < 2.25 

u i (t) = U 2 (t) = / .02 [l.3 sin (10 re t) + .4 sin (3 • 10 re t) + .25 sin (5 • 10 3* t)] 2.25 < t < 3.75 

(.32 [l.3 sin (l.9 re t)+ .4 sin (3 • 1.9re t) + .25 sin (5 • 1.9 re t)] i > 3.75 


This input is the Fourier series approximation of a square wave. Enough terms are retained 
such that the signal is sufficiently exciting. The frequency and amplitude of the square 
wave are also varied in order to test the robustness of the adaptation algorithm. 

Figures two and three compare the two system outputs with their corresponding 
estimates. The estimates appear to converge within one-half of a second. Changes in the 
signal have no apparent effect on the output estimate. 

Closed-loop simulations are obtained by using the MIMO minimum variance control 
law, and setting w(t) equal to the square wave given previously. Figures four and five 
compare the open-loop response (controller off), the closed-loop response, and the estimate 
of the closed-loop output. Clearly, the controller is effective in attenuating the time- varying 
noise disturbance. The closed-loop estimate does not converge on the actual system output 
because the exogenous input term, C(z)/A(z) w(t), is unmodeled by the controller. Since 
w(t) is unmeasurable, the coefficients of C(z) may not be determined using an algorithm 
such as RLMS. 



Figure two: Simulated Open-Loop Test of the RLMS Algorithm: 
A Comparison of Output One and Its Estimate 





Figure four: Closed-Loop Minimum Variance Controller Simulation: 

A Comparison of Controlled and Uncontrolled Output One and Its Estimate 









Open Loop Output 2 Estimated Output 2 Closed Loop Output 

Figure five: Closed-Loop Minimum Variance Controller Simulauon: 

A Comparison of Controlled and Uncontrolled Output Two and Its Estimate 

A MIMO adaptive minimum variance active noise controller has been implemented as a 
Turbo Pascal program on a 25 MHz 386-based PC. A Metrabyte DAS-20 data acquisition 
and control card provides the interface to the two microphones and speakers used in the 
experiment. A third speaker acts as a noise source, and is driven by a sinusoidal input. 
Noise cancellation experiments are conducted in a typical laboratory environment with the 
same speaker/microphone geometry shown in figure one. 

Figures six, seven, and eight illustrate the convergence of the Bn, Aj, and A 2 
polynomials from the RLMS algorithm, as well as the resulting error between the output 
and its estimate. It is important to note that the initial values for the coefficients have a 
very important role in the numerical stability and convergence of the estimation algorithm. 
This is not the case for the SISO controller. The MIMO RLMS performance is much less 
robust than its SISO counterpart, in that 22 coefficients must converge properly as opposed 
to the previous seven coefficients. This lack of robustness degrades both stability and 
performance. The change in controller effectiveness is obvious from the insertion loss 
plots given in figure nine. 

The MIMO controller exhibited the ability to adapt to changes in speaker/microphone 
geometry and acoustic impedance. Small changes in speaker and microphone placement 
did not affect performance significantly. Additionally, foam blocks were inserted in the 
acoustic path with no apparent change in the level of sound attenuation. 




Figure six: RLMS Adaptation of ARMAX Denominator Coefficients (A3 = 1) 



Figure seven: RLMS Adaptation of ARMAX Numerator 
Coefficients for Polynomial B 1 1 
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Figure eight: Error Between Measured Plant Output and Estimated 
Output from ARMAX Model 
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ENGINE MODEL DEVELOPMENT: 


A great deal of effort is spent by the research community in the area of gas turbine 
modelling. Typically, such models are highly complex, and represent the dynamics ot only 
one engine component, e.g. an axial flow compressor model [7]. Researchers attempt to 
describe system behavior with systems of partial differential equations, or more recently, 
differential equations which lead to chaos. 

While such models are useful for engine design and analysis, they are not appropriate 
for control design purposes. This need has motivated the development of state space 
models for control synthesis. One such model, HYTESS, is a four state model of a 
hypothetical twin-spool turbofan engine [8]. Proprietary state space models for specific 
engines also exist. Unfortunately, the simplicity of these models comes at a price - a lack 
of information concerning the internal pressure dynamics of the engine. The inclusion of 
internal pressures in the system state is crucial if the model is to be used for noise control. 
Clearly, a new state space model is needed. 

The bond graph technique pioneered by Paynter seems well suited for modelling 
turbomachinery, as it provides a uniform methodology for modelling mechanical, thermal, 
and fluid energy modes. Compressors and turbines are special cases of a general 
turbomachine in which interactions between the thermal and fluid domains result in shaft 
power (turbine), or compression of the fluid (compressor). Figure ten depicts this energy 
exchange as a simple bond graph [5]. 
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Figure ten: A Simple Turbomachine Bond Graph 


A typical gas turbine is composed of several components: a compressor, a combustor, a 
turbine, a shaft to transfer work from the turbine to the compressor, and an exit nozzle. In 
practice, turbofan engines are comprised of multiple compressor and turbine stages along 
with more than one driving shaft. The model developed here is for the simpler single stage 
case which appears as figure eleven. Extension to a more realistic engine representation is 
not difficult, as will be shown shortly. 




The energy flow, or power, associated with each component in a gas turbine engine is 
effected by one or more of the following mechanisms: 

• Resistance (Nozzles) - R 

• Capacitance (Fluid compressibility) - C 

• Inertance (Fluid flow) - 1 

Each component of the gas turbine engine may be described by a “basic functional unit” 
(BFU) bond graph comprised of these “R”, “C”, and “I” elements [6]. Since temperature 
and heat flow are used as energy variables, the resulting model is a pseudo bond graph. 
The model describes the relationships between pressure, mass flow, torque, and engine 
speeds as well as temperature and heat flow, and is shown in figure twelve. 



Figure twelve: Bond Graph Representation of a 
“Basic Functional Unit” 



The fundamental differential equations governing the generic BFU are found by 
applying constitutive laws between effort and flow variables. They are: 

£mi = f[Pti + DP t -Rmi] 

^- p to = ^-[m, - me - rho] 

£T to = ^-[rhi (T ti + DT ( ) - mo T lo - rile T lo + (Qor W)] 
at Cj 

The complete bond graph for the simplified gas turbine appears in figure thirteen, and is 
obtained by cascading multiple BFUs and eliminating unnecessary elements. 
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Figure thirteen: Bond Graph Representation of a Gas Turbine Engine 

The differential equations describing the complete bond graph are: 
d m _ 1 


^-ihci =p[P 0 - Rc nici] 


— P 

dt c 


_ 1 


C, 


•[n\i - (m tl - rhf)] 


cp 


^ [riici T 0 - (rhii - rhf) T co + q c T to] 
dt C c t 

^-Tbo = ^-[(m,* - rh f ) T co - rhti Tbo+ % AH C m f ] 

m,j = ^-[P co ■ Rr rh, j 

4r P.o = 7=r — [rhti - m to ] 
at L-tp 








4T to = ^ L [ni,iT bo - ihtoTto-ntTw] 

dt C, t 

^<o =t L[r-R 8 c] 

At first glance these eight state equations appear to be linear, however, the "R”, ‘T\ 
and “C” parameters are functions of the state as well as physical parameters associated with 
the gas turbine [6]. The parametric relationships are given as: 

j b L_ Component Length 

A Effective Component Cross-Sectional Area 

£ ^ LA _ Component Volume 

P n R T (Polytropic Index) (Gas Constant) (OuUet Temperature) 

£ m nCpLA P _ (Constant Terms) (Volume) (Outlet Pressure) 

* R (n- 1 ) y T (Constant Terms) (Outlet Temperature) 

Although the polytropic index, n, is lumped with other “constant terms,” it is important to 
note that it may also be a variable depending on the nature of the process. Namely, 
whether it is assumed isothermal or isentropic, n would change from a minimum to a 
maximum value. 

A linearized model of the system is obtained by selecting reasonable values for the 
physical parameters, selecting an operating point, and taking a Taylor series expansion 
about that point. The parameters used in the model are given in Table one. 


Parameter 

Value 

Parameter 

Value 

Lc 

! m 

nc 

0.8 

A c 

2.64 m 2 

m 

0.9 

Lt 

l m 

Tib 

0.98 

At 

2.64 m 2 

AHc 


L b 

0.5 m 

C P 

1005 J / k2K 

Y 

1.4 

n 

1.35 | 

R c 

UUliH 

Rt 


R s 

1.7 n^/s 

R 

287.1 J/ kaK 


Table one: Gas Turbine Parameters Used for 
Linear Model Development 




























Figure fourteen represents the responses of both the compressor and turbine outlet 
pressures to a step change in fuel flow. The pressure ratio of this engine appears to be 6: 1. 
Perhaps of greater interest, are the frequency response plots of the two pressures with 
respect to the fuel flow input shown as figure fifteen. A lightly damped pair of poles affect 
the response at 334 rad/sec in our single stage model. High noise levels in this frequency 
range have been reported in previous experimental studies. Simply adding additional 
BFU’s to the bond graph (representing additional turbine and/or compressor stages for 
example) could produce more peaks. 

The parameters used in the preparation of this model may not be representative values 
of actual engine parameters. The bond graph techniques shown here provide a convenient 
means of generating structured models of desired order, however, some means of 
parameter identification is needed to complete the model. Results from more sophisticated 
simulations or (preferably) actual engine frequency spectra could be used to complete this 
task. 



and Turbine Outlets tor a Step Change in Fuel Flow 
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and Turbine Outlet Pressures effected by the Fuel Flow 


CONCLUSIONS: 

A Multiple-Input/Multiple-Output adaptive minimum variance controller has been 
developed for active noise cancellation. The controller has been shown to be effective in 
simulation. The controller implementation was also successful at attenuating pure tone 
disturbances, and was able to adapt to changes in speaker/microphone geometry as well as 
acoustic impedance. The MIMO controller performance did not represent a significant 
improvement over the SISO controller in the experimental setup used during the 
implementation stage. MIMO controller performance appears to be limited by the lack of 
robustness of the RLMS algorithm used to identify system parameters. 

A framework for developing state space models of gas turbine engines suitable for 
acoustic controller design has been discussed. Bond graph model components are 
cascaded to represent multi-stage compressor and turbine systems. This modular approach 
allows a great deal of flexibility in determining appropriate model order and shaping the 
open loop characteristics. Some form of parameter identification is required before such a 
model would be of use. A hypothetical gas turbine model is presented in order to illustrate 
the types of results which could be obtained from a more sophisticated multi-stage model 
with properly identified parameters. 



SUGGESTIONS FOR FUTURE RESEARCH: 


The RLMS algorithm used for on-line parameter identification in the MIMO minimum 
variance controller exhibits a lack of robustness. This can lead to stability problems and 
poor attenuation. An alternate method of parameter identification using neural networks 
has been proposed by Chu and Shoureshi [9, 10]. Such techniques may improve both the 
stability and performance of an adaptive controller. 

The adaptive minimum variance control structure has been the focus of this study. An 
adaptive pole placement controller was also evaluated, but was abandoned due to the 
computational complexity it required. It is unclear if the conservative designs of robust 
control methodologies such as H^o and QFT are appropriate for noise control. 
Nevertheless, it would be of interest to evaluate these and other advanced control structures 
with regards to the noise cancellation problem. 

A multi-stage turbofan engine model is required for noise controller design and 
evaluation. The modelling approach presented here is meant to serve as a starting point for 
such research. Parameter identification and experimental verification of the model are 
important parts of this process. 

Finally, an important issue which must be addressed is the selection of actuation and 
sensing mechanisms for this application. Control inputs must be identified which satisfy 
the controllability requirements of the nonlinear, time-varying system of a jet engine. 
Similarly, the dual problem of sensor selection and observability must be considered. The 
harsh conditions which these sensors and actuators would be subjected to make this a 
challenging problem. 
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